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Abstract 

We calculate the full 0{as) radiative corrections to the three spin independent and 
five spin dependent structure functions that describe the angular decay distribution 
in the decay of a polarized top quark into a VF-boson (followed by the decay 

+ ui or by Q + q) and a bottom quark. The angular decay distribution is 

described in cascade fashion, i.e. the decay + Xf, is analyzed in the top 

rest system while the subsequent decay + ui (or Q + q) is analyzed 

in the W rest frame. Since the structure function ratios depend on the ratio mw/^t 
we advocate the use of such angular decay measurements for the determination of 
the top quark's mass. Our results for the eight 0{as) integrated structure functions 
are presented in analytical form keeping the mass of the bottom quark finite. In the 
limit mi, — > the structure function expressions reduce to rather compact forms. We 
also present results on the nifj = unpolarized and polarized 0{as) scalar structure 
functions relevant to the semiinclusive decay of a polarized top quark into a charged 
Higgs boson t(|) + X[, in the Two Higgs Doublet Model when mj, = . 



1 Introduction 



In the decay of an unpolarized or polarized top quark to the W-gauge boson and a bottom 
quark the is strongly polarized, or, phrased in a different language, the has a 
nontrivial spin density matrix. Furthermore, the spin density matrix of the W can be 
tuned by changing the polarization of the top quark. The polarization of the will 
reveal itself in the angular decay distribution of its subsequent decays + (or 

W+ ^q + q) Q. 

In the first stage one will aim to analyze the decay of unpolarized top quarks (or average 
over its polarization). The decay distribution of unpolarized top quark decay is governed 
by three structure functions which we shall refer to as Hu ( "unpolarized-transverse" ) , 
Hl ("longitudinal") and Hp ("forward-backward-asymmetric"). In fact, the CDF col- 
laboration has already presented some results on the measurement of the longitudinal 
component of the W based on the limited RUN I data The measurement has con- 
firmed the expected dominance of the longitudinal mode. The error on this measurement 
is quite large (~ 45 %) but is expected to be reduced significantly during RUN II at the 
TEVATRON to start in the spring of 2001. In RUN II one will produce (5 - 6) x 10^ top 
quark pairs per year and detector. This number will be boosted to 10^ — 10^ top quark 
pairs per year and detector at the LHC starting in 2006/2007. It is conceivable that the 
errors on the structure function measurements can be reduced to the 1 — 2 % level in the 
next few years 0. If such an accuracy can, in fact, be achieved and, having in mind that 
the 0{as) corrections to the top decay rate amount to 8.5 % 0, @, H, P, 0, ||, it is quite 
evident that one needs to improve on the known theoretical Born level predictions for 
the above three structure functions by calculating their next-to-leading order radiative 
corrections. 

At a later stage, when the data sample of polarized top quarks has become sufficiently 
large, one will be able to also analyze the decays of polarized top quarks. The top 
quark is very short-lived and therefore retains its full polarization content when it decays. 
Polarized top decay brings in five additional polarized structure functions which can be 
measured through an analysis of spin-momentum correlations between the polarization 
vector of the top quark and the momenta of its decay products. 

Polarized top quarks will become available at hadron colliders through single top 
production which occurs at the 33 % level of the top quark pair production rate 0. Future 

^ From this point on we shall drop explicit reference to the W'^ q q decay channel since it has 
the same angular decay distribution as — s- Z"*" + i^;. In fact the branching fraction into the two 
hadronic channels (d + w) and (s + c) exceeds that of the sum of the three leptonic channels by a factor of 
approximately two because of the colour enhancement factor. Although not explicitly mentioned further 
on, the existence of the hadronic decay mode of the is always implicitly assumed in the following. 
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e^e~ colliders will also be copious sources of polarized top quark pairs |TD|, |TT], |T^, 
p!4| , p!5| . For example, at the proposed TESLA collider one expects rates of (1 — 4) x 10^ 
top quark pairs per year. The polarization of these can be easily tuned through the 
availability of polarized beams (see e.g. [1^). Further, there is a high degree of correlation 



between the polarization of top and anti-top quarks produced in pairs either at e+e 
colliders [1^, 18, or at hadron colliders [^] which can be probed through the joint 



decay distributions of the top and the anti-top quark. 

In this paper we study momentum-momentum and spin-momentum correlations in 
the cascade decay process t — > + h followed by — > + 1^1. The step-one decay 
t W'^+h is analyzed in the t-rest frame where we study the spin- momentum correlation 
between the spin of the top and the momentum of the W . In step two we go to the rest 
frame of the W and analyze the correlation between the momentum of the lepton (or 
antiquark) and the initial momentum direction of the W . In technical terms this means 
we analyze the double density matrix of the decaying top quark and the produced W- 
gauge boson. This must be contrasted with the center of mass analysis of polarized top 
decay where the spin-momentum correlations are all analyzed in the rest system of the 
top quark (for an 0(0;^) analysis of this kind see |2^). Experimentally such a correlation 
measurement is easier, but from a theoretical point of view the cascade-type of analysis 
is advantageous because one can then better isolate the contribution of the longitudinal 
mode of the l^-gauge boson which is of relevance for the understanding of the electroweak 
symmetry breaking sector in the Standard Model. The results of the two analysis' are of 
course related through a Lorentz boost along the W direction. However, the azimuthal 
correlations to be discussed later are not affected by such a Lorentz boost and are thus 
identical in both types of analysis. 

The complete angular decay distribution is governed by altogether eight structure 
functions which we calculate analytically including their full 0(0;^) radiative corrections. 
One of the motivations for calculating the 0{as) radiative corrections is the fact that 
the radiative QCD corrections populate helicity configurations that are not accessible at 
the Born level. Take for example unpolarized top decay where, at the Born level, the 
cannot be right handed, i.e. cannot have positive helicity, due to angular momentum 
conservation when 771^ = 0. This implies that strictly forward production does not 
occur at the Born level. However, when radiative corrections are taken into account, 
right-handed W^s do occur and strictly forward production is allowed. As we shall see 
in Sec. 4 technically this means that the structure function combination {Hu + Hp)/ 2 
vanishes at the Born term level but becomes nonzero at 0{as) [0]. We shall, however, 
see that the 0{as) population of the right-handed W is rather small pS). The same 
statement holds true for the other structure function combinations that vanish at the 
Born term level. 
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In order to retain full control over the b mass dependence, and having also other ap- 
plications in mind, we have kept a finite mass value for the b quark in our calculation. 
This improves on our earlier calculation of polarized top decay where the b quark mass 
was neglected and where we limited our attention to the six (diagonal) structure func- 
tions that govern the polar angle distribution in the cascade decay The additional 
two (non-diagonal) structure functions calculated in this paper describe the azimuthal 
correlation of the plane of the top quark's polarization and the plane defined by the final 
leptons. In addition we determine the unpolarized and polarized scalar structure functions 
which are of relevance in the analysis of top decay into a bottom quark and a charged 
Higgs boson p3. We mention that our calculations have been done in the zero width 



approximation of the VT-boson. Finite width effects will be addressed in a forthcoming 
paper PH| (see also 



Most of the results in this paper are new. They have been checked against limiting 
cases and partial results obtained in other papers. We have checked our analytical 0{as) 
result for the total rate against the corresponding analytical rate result of Denner and 
Sack who also kept the b quark mass finite [Q. We find agreement. We took the zero 
6-quark mass limit of the six diagonal structure functions and obtained agreement with 
our previous results in |]16|. These had already been checked against the analytical results 
on the total rate obtained in ^ ^, |^ and on the longitudinal/transverse composition 
obtained in [^. All six (mass zero) diagonal structure functions had also been checked 
against the corresponding numerical results given in pTl ESI EO]. The unpolarized scalar 



structure function has been checked against the results of p4| 



The central topic of this paper is the analysis of polarized top decay. We therefore 
mostly limit our attention to results valid in the limit rrib — > in the main part of our 
paper. This leads to enormous simplifications in the analytical rate formulas. The quality 
of the rrih = approximation may be judged from the Born term rate which increases by 
0.27% going from nib = 4.8 GeV to = 0. The full nib ^ structure is given in Sec. 8 
and the Appendices. Apart from retaining full control over 7^ effects the finite mass 
results are needed e.g. in the theoretical analysis of semileptonic 6 — >■ c decays where the 
c-quark mass can certainly not be neglected. 

Our paper is structured as follows. In Sec. 2 we define a set of three spin independent 
and five spin dependent structure functions through the covariant expansion of the decay 
tensor resulting from the product of the two relevant current matrix elements. The 
eight invariant structure functions are related to eight helicity structure functions which 
form the angular coefficients of the angular decay distribution. In order to facilitate 
the calculation of the tree graph contributions we define a set of five covariant projection 
operators and a covariant representation of the spin vector of the top. These projectors can 
be used to covariantly project the requisite helicity structure functions from the hadron 
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tensor. The advantage is that one thereby obtains the appropriate hehcity structure 
functions and scalarizcs the tensor integrands needed for the tree graph integration in 
one go. In Sec. 3 we derive the exphcit form of the angular decay distribution in terms 
of the eight hehcity structure functions for top decay. We also specify the changes in 
the angular decay distribution needed for antitop decay. Sec. 4 contains our Born term 
results. In Sec. 5 we list our results for the nib — one-loop contributions. In Sec. 6 we 
provide expressions for the 0(0;^) tree graph contributions and discuss technical details of 
how we have handled the necessary tree graph integrations. We mention that the infrared 
divergencies are regularized by a finite small gluon mass. In Sec. 7 we take the mf, — > 
limit of the 7^ results in Sec. 8 and present rather compact analytical 0{as) formulas 
for the various structure functions. Sec. 7 also contains our numerical results in the mi, — 
approximation. Sec. 8 gives our analytical results on the tree graph integrations plus the 
one-loop contributions for 7^ 0. Sec. 9 provides a summary and our conclusions. In 
particular, we emphasize that angular measurements as advocated in this paper can be 
utilized to measure the mass of the top quark. In Appendix A we provide a complete list of 
m{, 7^ basis integrals that appear in the calculation of the tree graph contributions. This 
set of basis integrals should also be useful for other 0{as) or 0{a) radiative correction 
calculations. The requisite coefficient functions that multiply the basic integrals in the 
structure function expressions are listed in Appendix B. Appendix C, finally, contains the 
one-loop contribution in the 7^ case. 

2 Invariant and helicity structure functions 

The dynamics of the current-induced i — > 6 transition is embodied in the hadron tensor 
H'^^ which is defined by 

Xb 

X {t{pt,st)\r+\Xb){Xb\J^\t{pt,st)) , (1) 

where dUf stands for the Lorentz- invariant phase space factor. In the Standard Model 
the weak current is given by — qhl^Ph^t with Pi = \{1 — 75). 

We are working in the narrow resonance approximation of the W^-boson and set q^ = 
Tfiy\^ as indicated in the argument of the hadron tensor. Thus the hadron tensor is a 
function of the energy go of the W alone. Since we are not summing over the top quark 
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spin the hadron tensor also depends on the top spin St as indicated in Eq.(|T]). The 
structure of the hadron tensor can be represented by a standard set of invariant structure 
functions defined by the expansion 

H'^^={- g'^^ H, + pfp,^ H2 - le^^'^t^pq. H,) + 

- (g-s,)(-(7^^Gi+prprG'2-^e^^^>pg.G3) + (2) 

where the (i=l,2,3) and Gi (1=1,2,3,6,8,9) denote unpolarized and polarized structure 
functions, respectively. 

In the expansion (Q) we have kept only those structure functions that contribute in 
the zero lepton mass case. We have thus omitted covariants built from and/or q". 
We have also dropped contributions from invariants that are fed by T-odd or imaginary 
contributions which are both absent in the present case. 

In the expansion one has still overcounted by one term since there is a relationship 
between the three parity conserving (p.c.) spin dependent covariants appearing in (H) due 
to the identity of Schouten. The identity between the three covariants reads 

q-St e^^^^.pg. - g'e^^^^.pSi, + q-pt e^'^^q^St,, = (3) 

We shall, however, keep the overcounted set of nine invariant structure functions in (0) 
for reasons of computational convenience. 

In this paper we shall only be concerned with two types of intermediate states in (|l]), 
namely \Xh) = \b) (Born term and 0{as) one-loop contributions) and \Xh) = 16 + 5') 
{0{as) tree graph contribution). The Feynman diagrams contributing to the respective 
processes are drawn in Fig. 1. 

The angular decay distribution that we are aiming for is given in terms of a set of 
angular decay coefficients which are linearly related to the set of unpolarized structure 
functions and polarized structure functions Gi defined in Eq. (|]). The relevant linear 
combinations are given by 



Hu = H++ + H__ = Hu + H22, (4.a) 

Hl = Hoo = H^3, (4.b) 

Hp = H++-H^^=t{Hi2-H2i), (4.c) 

Hup = H^44) + H_4s\)=Hnis\)+H22i4), (4.d) 
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Hlp = Hoo{s\) = Hs3{s\), (4.e) 
Hp. = H^M)-H-{s\)=^{H,,{s\)-H^M)). (4.f) 



2V2 
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^ (H.sisT) + H.^isf)), (4.h) 



Ha^ = ](hMsT) + Ho4sT) + HMsT) + Ho4sT)) (4.i) 



2V2 



{H,,{sr) - Hs2{sm (4.j: 



where Hx^.y^ = H^ue*^{\w)^^{^'w) ^^e helicity projections of the polarized and 
unpolarized pieces of the structure functions H^'^. The e*^{Xw) and e^{\w) are the usual 
spherical components of the polarization vector of the W gauge boson. In the top quark 
rest system with = (go; 0, 0, |g |) and |g | = {qq — m^)^/^ they read 

e'^(0) = ^(|g|;0,0,go), (5) 
rriw 

e'^(±) = -^(0;Tl,-^,0). (6) 



In Eqs.( [4.a| -^7][) we have also included the Cartesian components of the helicity structure 
functions in the ly-boson rest frame which are useful to have for some applications. For 
notational convenience we shall often refer to the set of helicity structure functions by 
their generic names. Thus we shall frequently use U for Hu and for H^p, etc.. 

The rest frame components of the longitudinal ("1") and transverse ("tr") polariza- 
tion vector of the top are simply given by s{ = (0; 0, 0, 1) and s*^ = (0; 1, 0, 0). For the 
unpolarized helicity structure functions one sums over the the two diagonal spin con- 
figurations of the top while one takes the differences of these for the polarized helicity 
structure functions (in the 2-basis for s' and in the x-basis for sf'). When computing the 
polarized structure functions from the relevant Dirac trace expressions one thus has to 
replace {-^t + ^t) in the unpolarized Dirac string by {"^t + ^t)(l + 75;^*)- Note that the 
longitudinal component contributes only to the diagonal helicity structure functions U, L 
and F while the transverse component contributes only to the non-diagonal structure 
functions I and A. The physics behind this will become clear when we write down the 
angular decay distribution in Sec. 3. 

It turns out that it is rather convenient from the computational point of view to 
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represent the helicity projections in (^7a|-^) (defined by the gauge boson polarization 
vectors and the top polarization vector) in covariant form. One has 



=i/^.Pr z = U,L,F, (7.a) 

H,P=H,,{4)Fr i = U,L,F, (7.b) 

= H,,{sf)]Pr 1 = 1, A. (7.c) 

The covariant projectors onto the diagonal density matrix elements are given by 

Pt-q u\( u Pt-q 









= -j"" + - 


P^" 




mt \q \ 



[Pt - ^^q^) {pi - -::j-q^ 



m 



2 ' 
W 



(8.b) 



Pt,aqf3, 



where e°^^^ = —1. We do not write out the projector for the unpolarized-transverse 
component U but note that it can be obtained from the combination P(/+i — F'^ . 

The projectors onto the transverse-longitudinal non-diagonal density matrix elements 
are given by 

2v2 mt \q\^ ^ ' > 

P7 = --^^^{.e'^"^-6„(xK,g,(pr - ^g^) - /i - z^}. (S.e) 
2v2 m^t \qY ^ ^ vtii^- ' > 



w 



They involve the the transverse polarization vector of the IV-gauge boson ej^x) = (0; 1, 0, 0) 
pointing in the x-direction. 

The covariant representation of the longitudinal component of the polarization vector 
of the top spin vector s\ is given by 

sY = Uq'-'^P^, (9) 
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whereas its transverse component s*^ reads 

sr''^ = (0; 1,0,0). (10) 



Note the inverse powers of \q\ = \Jql — that enter the L, T, F, / and A projectors 
and the longitudinal polarization vector. They come in for normalization reasons. These 
inverse powers of |g | will make the necessary tree graph integrations to be dealt with in 
Sec. 6 and in the Appendices A and B somewhat more complicated than the total [U + L) 
rate integration which has a rather simple projector as Eq. (8.2) shows. 

As mentioned in the Introduction, the covariant forms of the projection operators ( ^.a| - 
P^ ) and the polarization vectors (^ and ( pUj ) are quite convenient for the calculation of 
the 0{as) tree graph contributions to be dealt with in Sec. 6 . The covariant projectors 
allow one to scalarize the tree graph tensor integrands and to project onto the requisite 
helicity structure functions in one go. 

Although we shall mostly work in the helicity representation of the structure functions, 
it is sometimes convenient to have available the set of linear relations between the helicity 
and invariant structure functions. These can easily be worked out from the expansion (H), 



the projectors (p.a|-g.e) and the polarization vectors ([101). One has 



Hu 


= 2ifi. 


(ll.a) 




= m^Hi + \q\^mlH2, 


(ll.b) 


Hf 


= 2\q\mtH3, 


(ll.c) 


Hup 


= 2|g|Gi, 


(ll.d) 




= |g|(m^Gi + Igpm^Ga - 2qomtGQ), 


(ll.e) 


Hpp 


= 2 \q\'^mt G^-2mt Gs -2qoGQ, 


(ll.f) 


Hjp 


V2 mw 


(11-g) 


Hap 


1 nitqo ^ 1 ^ 
V2 mw V2 


(ll.h) 



Note that the three structure functions G^, Gg and Gg always contribute in the two 
combinations (m^Gs + Gg) and (go^iGa — Gg) proving again that there are only eight 
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independent combinations of structure functions. If desired, Eqs.( |ll.a| - p.l.i:]| ) can be in- 
verted such that the invariant structure functions can be expressed in terms of the hehcity 
structure functions. The inversion has to be done in terms of the two above hnear com- 
binations of G3, Gs and Gg. Since our later results will always be presented in terms of 
the helicity structure functions, we shall not write down the inverse relations here. 



3 Angular decay distribution 

We are now in the position to write down the full angular decay distribution of polarized 
top decay into and b followed by the decay of the the into (Z^ + ui). As noted 
before, the full angular decay distribution of the decay t(t) — ^ W^{—^ + ui) + X^, 
including polarization effects of the top quark, is completely determined by the three 
unpolarized and the five polarized helicity structure functions. Although the necessary 
manipulations to obtain the angular decay distribution involving Wigner's D;^ 



mm' 



functions are standard (see e.g. |]30|), it is quite instructive to reproduce the results here. 
To this end, it is useful to define helicity structure functions H^'^^y where the helicity 
label of the top quark is made explicit. Put in a different language the four-index object 

At A' 

^Xw^' unnormalized double density matrix of the top and the W . The double 

density matrix is Hermitian, i.e. it satisfies 



= (^It)- (12) 



As has been remarked on before the elements of the double density matrix are real in 
the present application. The double density matrix is therefore symmetric. The relation 
of the components of the double density matrix to the previously defined unpolarized and 
polarized helicity structure functions is given by 



(13.a) 
(13.b) 
(13.C) 
(13.d) 
(13.e) 



Hu = 


-"++ 


+ H-- + H++ 


+ HZZ 


Hl = 


^00 






Hf = 


Htt 




-HZZ 


Hup = 


nx 




-HZZ 


Hlp = 


^00 


~ ^00 ' 





10 



Hpp = ~ -^++ ~ H-~ + H_ 



1 





(13.f) 


+0 H_q), 


(13.g) 




(13.h) 



Hip = ^iH_^_Q + H^^ - H_q - Hq_ ) - -, 
Hap = \iH+^ + H^+ + HZ^ + H+Z) = \\ 

For ease of notation we have used (±)-labels for both the hehcities of the top (At = 
±1/2) and the transverse hehcities of the W gauge boson {\w = ±1). In the case of 
the non-diagonal structure functions Hjp and H^p one can make use of the fact that 
the double density matrix is symmetric (for real coefficients !) to simplify the structure 
functions as indicated in the last two lines of Eqs. (|13.g| - [T3li| ). From the fact that we 



are not observing the spin of the Xf, system in our semi-inclusive measurement one has 
Ajsfi, = A'^^ leading to the constraint \w — A'^ = \t — A'^. From this constraint it is 
immediately clear that the polarized structure functions L and F are associated with 
the longitudinal spin of the top and the structure functions / and A are associated with 
the transverse spin of the top. 

The angular decay distribution can be obtained from the master formula 



w 



iy(^P,^,0) cx Y: e^^'""''-^^4^iW4wiW^AtA' Pa.a^(M, (14) 



-''w — A(^— At — A 

where P\t\[{Op) is the density matrix of the top quark which reads 



, . 1 fl + PcosOp PsinOp \ 
PA.A^(^p)- 2 pgin^^ 1-Pcos^pj- ^^^^ 



P is the magnitude of the polarization of the top quark. The sum in Eq. (jM]) extends over 
all values of Xw, A^y, Xt and A^ compatible with the constraint Xw — X'^r = Xt — X[. The 
second lower index in the small Wigner (i(^)-function c^i^i is fixed at m = 1 for zero mass 
leptons because the total m-quantum number of the lepton pair along the direction is 
m = 1. Because there exist different conventions for Wigner's rf-functions we explicate 
the requisite components that enter Eq. (p^ : d\i = (1 + cos6')/2, dl^ = sm9/\/2 and 
di^-^ = (l-cos^)/2. 

Including the appropriate normalization factor the four-fold decay distribution is given 

by 



1~1 a~l ^ = ^ 7^ g {- i/£/ + Pcos^pi/[;p 1 + cos 6) 

dqod cos Opd cos Odcp 47r v2 7r '^8 
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+ ^{Hl + PcosOpHLp) sin^e + ^{Hf + PcosepHpp) cose (16) 

3 3 

H -^P sin epHjp sin 26' cos (j) H — -^P sin epHj<^p sin cos > 

2\/2 \/2 ^ 

We take the freedom to normalize the differential rate such that one obtains the total 
t h + rate upon integration and not the total rate multiplied by the branching 
ratio of the respective decay channel. 

The polar angles and ^, and the azimuthal angle that arise in the full cascade- 
type description of the two-stage decay process t{]) — >• W^{-^ + ui) are defined in 
Fig. 2. For better visibility we have oriented the lepton plane with a negative azimuthal 
angle relative to the hadron plane. For the hadronic decays of the W into a pair of light 
quarks one has to replace {l^^vi) by (g, g) in Fig. 2. We mention that we have checked 
the signs of the angular decay distribution Eq. (|1^) using covariant techniques. 

As Eq. ([T6|) shows the non-diagonal structure functions Hip and Hj^p are associated 
with azimuthal measurements. This necessitates the definition of a hadron plane which is 
only possible through the availability of the x-component of the polarization vector of the 
top (see Fig. 2). This is the physical explanation of why the two structure functions Hip 
and Hj^p are functions only of the transverse component of the polarization vector of the 
top quark. For similar reasons the polarization dependent structure functions Hup^H^p 
and Hpp depend only on the longitudinal component of the polarization vector. 

Setting P = in Eq. (^^ one obtains the decay distribution for unpolarized top decay. 
If desired, the transverse part of the unpolarized angular decay distribution can also be 
sorted in terms of decays into transverse-plus and transverse-minus VT-bosons given by 
the structure function combinations [U + F)/2 and [U — F)/2 which multiply the angular 
factors (1 + cos 6')^ and (1 — cos^^)^, resp., as done e.g. in p3[] . 

If there were an imaginary part in the one-loop contribution one would have two 
additional contributions to the angular decay distribution proportional to sin (p. This can 
be easily seen with the help of Eq. (|T^. We concentrate on those terms in the angular 
decay distribution that are proportional to the off-diagonal terms p_|__ in the density 
matrix of the top. The relevant terms read 

H+o e+''^ + Hq + e-''*' = 2 {Re{HXo) cos0 - Im(i7+o ) sin0), (17.a) 
^o- e'^'^ + HZo e'"'*' = 2 {Re{H+~) cos (j) - lm{H+:) sin 0) . (17.b) 
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The real contributions multiplying the angular factor cos have been included in the 
angular decay distribution ([T6|) while the imaginary part contributions lm{H^Q) and 
lm{HQZ) multiplying sin0 do not appear in Eq. (|T6|) since the 0{as) contributions cal- 
culated in this paper are purely real. The helicity structure functions lm{H^Q) and 
Im(ifQr) are conventionally called T-odd structure functions and are contributed to by 
the imaginary parts of loop contributions and/or by CP-violating contributions which, as 
has been emphasized before, are not present in this calculation. 

Of interest is also the corresponding angular decay distribution for polarized anti-top 
decay W~{-^ l~ + vi) + X5. The angular decay distribution is changed due to the 

fact that the total m-quantum number of the lepton pair in the /~ direction is now m = — 1. 
The relevant components of the small Wigner (i-function are now d\_i = (1 — cos9)/2, 
dl_i = —sm9/^/2 and d\^_^ = (1 + cos6')/2. This can be seen to result in a sign change 
for the angular factors multiplying the F, and terms (and no sign change for the 
other terms). The structure functions of anti-top decay are related to those of top decay 
by CP- invariance. The p. v. structure functions F, U^, and will undergo a sign 
change whereas the p.c. structure functions U, L, F^ and keep their signs. Overall 
this means that the signs of the unpolarized terms in Eq. (|16D will not change their signs 
while the polarized terms will change signs when going from top decay to anti-top decay. 
To be quite explicit, if one wants to use the results of this paper to describe anti-top 
decay, the only required effective change is to change the signs of the terms multiplying 
the U^, L^, F^, and A^ structure functions in the angular decay distribution Eq. (|16|), 
using, however, the same structure functions as written down in this paper. 



4 Born term results 

The Born term tensor is calculated from the square of the Born term amplitude (see 
Fig. 2(a)) given by 

M^ = \4,-^M,7'^^(l-75)«*. (18) 

We omit the coupling factor Vtbg/V2 = 2mwVtb{G p / V2y^'^ and write for the Born term 
tensor (the spin of the b quark is summed) 

= ^Tri^t + mk)Y{l - I5M + mt){l + l5/st)Y{l - 75). (19) 
Since only even- numbered 7-matrix strings survive between the two (1 — 75)-factors 



in ([T9[) one can compactly write 



13 



B^^ = 2(prPr + P'tPl - g'^Pt-Pb + le^^'^^PbM), (20) 

where 

Pt =Pt - mst- (21) 

It is not difficult to obtain the Born term hehcity structure functions from (^^. This 
can be done in two ways. One can either read off the invariant structure functions ac- 
cording to the covariant expansion Eq. (^. The nonvanishing elements are given by 
= m^(l — + y^), Bh2 = —'^Bhs = 4 for the unpolarized invariants and by 
Bg^ = Bqq = -Bgs = ~Bgc) = — 2mt for the polarized invariants (the notation is self- 
explanatory). These can then be converted to the helicity structure functions using the 
linear relations ( |ll.a| - pni|) . Or, one can directly compute the helicity structure functions 
from (|20| ) by using the covariant projectors defined in Sec. 2 (c/. Eq.(p^ p^)). 

In order to find the relation of the Born term tensor B^^'^ to the hadron tensor H^^ 
defined in Sec. 2 one has to insert the appropriate one-particle 6-quark state into Eq. (|l|) 
and then one has to do the requisite one-particle phase space integration. Technically this 
is done by rewriting the one-particle phase space as 

jdlib = = jd%5{pl-ml). (22) 

One can easily do the four-dimensional rf^ph integration in Eq.([^) with the help of the 
four- dimensional (5-function 5^{pt — q — Pb)- This converts pi in the argument of the 5- 
function in Eq.(p2|) into {pt — qY- Rewriting the argument of the 5-function in terms of 
go one finally arrives at 

Ami 2 rrit 

We will present our results in table form where we use the scaled variables x = niy/ /irtt 
and y = mi,/mt as well as the abbreviation \q \ = (mt/2)\/X with A = A(l,x^,?/^) = 
1 + + — Ix^y^ — 2x^ — . The first column in Table 1 contains the ^ 0, or 
equivalently, y 7^ results. In the second column we have set mf, = (?/ = 0). In order 
to assess the quality of the m^, = approximation for the various rate functions we have 
listed the percentage increments when going from the mi, 7^ case to the mb = case 
including the phase space factor |g | that multiplies the helicity structure functions in the 
rate formula Eq. (0). In this comparison we have used mb = 4.8 GeV together with 
mt = 175 GeV and = 80.419 GeV. The increment due to the phase space factor 
|g| alone amounts to 0.15%. Note that one may have overestimated the mass effect 
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since a fixed pole mass, rather than a running mass which is smaller at the high scale of 
the top mass, is used. For example, taking one-loop running and the same bottom pole 
mass as above one has fhb{mt) = 1.79 GeV. The increment in the total rate going from 
fTT'b{mt) = 1.79 GeV to m;, = would then only be 0.04 % as compared to the 0.26 % given 
in Table 1. 

In the rrih = case listed in column 3 of Table 1 one observes the simple patterns 
Bjj = —Bijp = —Bp = Bpp, Bp = Bpp and Bjp = —Bj^p. This pattern results from 
the fact that a massless 6-quark emerging from a {V — A) vertex is purely left-handed. 
Since from angular momentum conservation one has At = Xw — Xb with Xb = —1/2 one 
has the constraint A* — Xw = 1/2. This implies that only the helicity configurations 
{Xt = —l/2;Xw = —1) and (A^ = +1/2; Xw = 0) are non-vanishing. A quick look at 
the relations ( |13.a| - P^XTiD allows one to readily verify the rrib = pattern in Table 1. 
For nib ^ there is a leakage into right-handed bottom mesons resulting in a breaking 
of the above pattern as can be observed in the rrib ^ column of Table 1. As noted 
in the Introduction these simple patterns are also not valid at 0{a^ even for massless 
bottom mesons because of the additional gluon emission including an anomalous spin-flip 
contribution [^. When the relevant nib = Born term helicity structure functions from 
Table 1 are substituted in (|16D we reproduce the angular decay distribution as written 
down in [|I7| . 

For completeness we have also included the two Born term scalar helicity structure 
functions Bs and Bgp in Table 1. They are obtained by use of the scalar projector 
IPs' = <l^(f /fT^'w- That they are identical to their longitudinal counterparts Bp and Bi^p 
even for 7^ is a dynamical accident specific to the Born term level and does not 
hold true in general as e.g. evidenced by the 0{as) contributions to be discussed later 
on. These become equal to each other only in the limit — > 00 as will be discussed in 
Sec. 7. The ntb 7^ Born term equalities Bp = Bijp and Bu = Bpp can be seen to result 
from the fact that the double density matrix elements H^l^ and H^t vanish at the Born 
term level due to angular momentum conservation (see p.3.a| ). 

In Fig. 3 we present a lego plot of the two-fold {mb = 0) Born term angular decay 
distribution in cos 6 and cos 6p which results after taking the azimuthal average of Eq. (|T6|) . 
We have divided out the total Born term rate from the differential rate resulting in the 
hatted differential rate distribution as defined in Eq.(|5^). We have set P = 1 in Fig.3. 
The lego plot shows that the cos^ and cosOp variation of the two-fold angular decay 
distribution around its average value of 0.25 is quite strong. This will facilitate the 
experimental measurement of the structure functions Tu, Tp, Tp, Tup, Tpp and Tpp . 

Finally, for the sake of definiteness we list the Born term rate in terms of the Born 
term function Bjj^p. One has 
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Born 

term 


m,, ^ 


mjy = 


increment 




,1^2 ^ { { -\ ^,2\2 1 ^2/1 o ^2 1 ^ ,2\ \ 

m^ — [[l-y ) +x [l-Zx +y )) 
x-^ 


,1^2 / 1 ^2\ / -I 1 o^2\ 

m^ — (1— a; j(l + 2a; ) 

X 


+U.z7 70 


I ) 


"^t V A— (l-2x -y j 


™2 -'^ /I ^2\/-i <-)^2\ 

X 


+\JAZ 70 


Bu 


0™™2/-i ^2 1 ^,2^ 

2mt{l-x +y ) 


0™™2/'-i ™2\ 

zm^ [1 — x ) 


+(J.(J5 70 




— zm^ V A 


o,™. 2 / 1 ™2\ 

— zm^ (1— x j 


1 n on 0/ 

+(J.29 70 




m^—[[l-y ) -X (1+2/ )) 

X 


^2 1 /I 2\ 


+U.OD 70 




X 


^2 -'- /I 2\ 


+U.d7 70 


Bp 


—Zm^ V A 


-2m^{l-x ) 


1 n on 0/ 
+U.zy /o 


Bpp 


2mt(l-a;^+7/^) 


2mi^(l-j;^) 


+0.05 % 


Bs 




m?l(l-x^) 


+0.36 % 


Bsp 


m,2^/A4(l-2/') 

X 




+0.37 % 


Bjp 


1 1 

— V2m^VX- 


-lV2m?l(l-x^) 


+0.29 % 


Bap 


iv^m?l(l-x^-y^) 

^ X 


iy2m?i(l-x^) 

Zi X 


+0.24 % 



Table 1: Born term helicity structure functions Bi {i = U+L, U^+L^, U, U^, L, L^, F, 
F^, S, S^, I^, for rrih ^ and = 0. Fourth column gives the percentage increment 
when going from 7^ to = including the phase space factor 
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G 



pm^ \q\ 



AV2 



vrmr 



U+L- 



(24) 



5 One-loop contribution 

The one-loop contributions to fermionic {V — A) transitions have a long history. Since 
QED and QCD have the same structure at the one-loop level the history even dates back 
to QED times. 

Our reference will be the work of Gounaris and Paschalis plf (see also |^) who used 
a gluon mass regulator to regularize the gluon IR singularity. The one-loop amplitudes 
are defined by the covariant expansion (Jj = qblfilt, = Qblfil^Qt) 



{b{pt)\J]^\t{pt)) = uM{i^F^ +Pt,^F^ +Pb,^F^}ut{pt), (25.a) 
{KPb)\J^\t{Pt)) = Ub{Pb){l^F,^+Pt,,F,^+Pb,^F,^}lMPt)- (25.b) 



In the Standard Model the appropriate current combination is given by — J^- 

We shall immediately take the limit — >^ of the one-loop expressions given in 
(see also Appendix C)0. Keeping only the finite terms and the relevant mass (M) (Iny 
and \yi' y) and infrared (IR) (In(A^)) singular logarithmic terms one obtains the rather 
simple result 



''." = l-^C'K*+^Ml--^)-Hn(^^)+ (26.a) 

+2'°(fT^)'"(Tr;;i)+2Li.(.^)), 

FV = -F- = ^^C.4( + 1 + i^lnd - (26.b) 

rrit 47r ^ ' 



^ We have recalculated the one-loop results of Ref. |34| and have found an acknowledged typo in the 
scalar form factors J3((3^) and H^jQ^) of The typo is corrected by replacing the factor (m2 — mi)/Q'^ 
in the last line of Eq.(A.8) of Ref. |||| by (ma - mi)/{2Q'^). 
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where we have denoted the scaled gluon mass by A = mg/rrit. The dilog function Li2(x) 
is defined by 



U,{x):^- l^^^^^-^dz (27) 
J z 



Note that the one-loop contribution is purely real. This can be understood from 
an inspection of the one- loop Feynman diagram Fig. 1(b) which does not admit any 
nonvanishing physical two-particle cut. The fact that one has — Ff^ and F-^ = —F/^ 
for i = 2,3 results from setting the 6-quark mass to zero. This can be seen by moving 
the chiral (1 — 75) factor in the one-loop integrand numerator to the left. Because nib 
is set to zero the Dirac numerator string will thus begin with ^^(l + 75) leading to the 
above pattern of relations between the loop amplitudes. We mention that the gluon 
mass regulator scheme can be converted to the dimensional reduction scheme by the 
replacement log — > 1/e— 7E+log47r/i^/g^ where 2e = A — N, is the Euler-Mascharoni 
constant 7^ = 0.577 . . . , and ji is the QCD scale parameter. 

6 Tree graph contribution 

The tree graph contribution results from the square of the real gluon emission graphs 
shown in Fig.lc and Id. Omitting again the weak coupling factor Vtbg/\/2 for the time 
being the corresponding hadron tensor is given by 



{PfPt){k^pl + k^pt - k-pbg'"' - i e'^^'^'^k^pb,/,) + 



+ 



^ k-Pb 



+ 



k-pt 

- {Pfk){{pt-kYpl + {pt-kYp^b - {Pt-k)-Pbg^'' -ie''^^''{pt-kYPb,p) 

- iPt -Pb) {k^ Pl + k'fb-k-Pb g'" - i e'^^^'^ka p^) + (Pt -Pb) {k' P'^ + k'' pf - k-Pt g^") + 

- {k-Pb){ptPt+PtPt-PfPtg'"') + {k-Pt){{Pb+kY p\+{pb+kY pt+{pb+k)-Ptg^''') + 
+ {k-Pt) {2p^,pl - Pb-Pbg^"") lif^^"" (k-Pt) + e^P-^^ k-'Pt,, - e^^^^ k^Ptn)Pb,aPt,0 + 
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>7 



■SGF 



A 



Anas Cp 







Pb-Pt 



) 



(29) 



SGF ■— — 



(k-pb) 



(k-pt) 



{k-pb){k-pt) 



where k is the 4-momentum of the emitted gluon. Asgf is the IR-divergent soft gluon 
function and pt = Pt — nitSt as in Sec. 4 0. 

We have isolated the IR-singular part of the tree-graph contribution by sphtting off 
a universal soft gluon factor which multiplies the lowest order Born term tensor B^^ . 
This facilitates the treatment of the soft gluon singularity to be regularized by a (small) 
gluon mass mg. Since the soft gluon factor is universal in that it multiplies the lowest 
order Born contribution, the requisite soft gluon integration has to be done only once 
and is identical for all eight structure functions. The result for the integrated soft gluon 
function is given in Sec. 8. Integrating only the soft gluon function /S.sgf and neglecting 
the finite part in Eq.(|28|) amounts to what is called the soft gluon approximation. We 
emphasize that we always include the full tree-graph contribution (soft plus finite part) 
in our calculation. Also, we integrate over the full phase space of the gluon, and not only 
up to a given energy cut-off of the gluon. 

We have deliberately used a calligraphic notation for the tree graph hadron tensor 
Ti.^'^ in (|28| ) since Ti^'^ is not the hadron tensor H^^ defined in Sec. 2. In fact, the mass 
dimension of Ti^^ differs from that of H^^ . To relate the two hadron tensors one has to 
do the appropiate phase space integration on the tree graph hadron tensor. 

Next one makes use of the covariant projection operators and the covariant forms 
of the longitudinal and transverse polarization vectors defined in Sec. 2 to obtain the 
contributions to the three unpolarized and five polarized structure functions. Since we are 
aiming for a fully inclusive measurement regarding the system the resulting expressions 
have to be integrated over the full two-dimensional phase space. As phase space variables 
we take the gluon energy and the W energy go where the k^ integration is done first. 
The phase space limits of the respective integrations are given by 



■^Contrary to the Born term case the polarization of the top quark cannot be accounted for by replacing 
all Pt momenta by their barred counterparts pt ■ 



ko- <ko< A;o,+ 



(30) 



and 



mw <qo< 



2mt 



(31) 
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where 



{m - qo){Ml - 2qQmt) ± \l ql - mly^{Ml - 2qQmtY - 'iml'inl 

^o,± = 7T? — ~, 2 n \ (^^) 

2{mt + - 2qoTnt) 

and 

Ml := + ml^, - ml ± ml. (33) 

It is clear from Eg. ( ^0] - ^ ) that the integration boundaries considerably simplify when 
the gluon mass is set to zero. In particular the second square root factor in the fco,± 
boundary turns into a polynomial in qo which is an essential simplification for the second 
go integration. This observation is at the core of our tree-level integration strategy exem- 
plified by the partitioned form of Eq. ([2^). The soft gluon singularity has been isolated 
and brought into a simple form. The remaining part of the tree-graph contribution is IR 
finite and can be integrated without the gluon mass regulator. 

The integration over the gluon energy ko (fco - < < ^o,+) is simple and the results 
will not be presented here in explicit analytical form. Instead we present some represen- 
tative results on the differential W^-boson energy distribution that result from the real 
gluon emission graphs Fig. Ic and Id in graphical form in Figs. 4 and 5. Fig. 4 shows 
the VT-boson energy distribution for the total rate dVu+hl dq^. The energy distribution 
rises sharply from the lower energy limit, where the M^-boson is produced at rest, then 
increases rapidly over the intermediate range of VT-boson energies and finally rises sharply 
again towards the end of the spectrum, where the soft gluon singularity is located. In 
Fig. 5 we show the same distribution for the partial rate into positive helicity W-bosons 
dT^/dqo (r_|_ = \{Tu + Tp)) for = and for ^ 0. As mentioned before there 
is no Born term contribution to dT^/dqo for mj, = and thus dV^/dqa possesses no IR 
singularity in this limit. The absence of the IR singularity in the mt = case (dashed 
line) is quite apparent in Fig. 4. The distribution rises moderately fast from the lower end 
of the spectrum, then turns down over the intermediate range of energies and finally tends 
to zero at the end of the spectrum where the phase space closes. The mh = (dashed 
line) and m^ ^ (full line) distributions lie on top of each other for most of the lower 
part of the spectrum. Starting at around 4.8 GeV below the upper phase space boundary 
the two distributions begin to diverge from each other. Whereas the = curve turns 
down and goes to zero at the end of the spectrum the mb curve starts to rise again 
and, in fact, tends to infinity at the end of the spectrum due to its IR singular behaviour. 
Note the huge differences in scale of the dTu^L/dqo and the dT^/dqo distributions which 
will be reflected in big differences in the total a^-corrections for the two respective rates. 

The second integration over the energy of the ly-boson is more difficult. Details can 
be found in Sec. 8 and in the Appendices. As it turns out the analytical m^ ^ results are 
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quite lengthy. We thus chose to present our rrih = results first since they are sufficiently 
simple to be presented in compact form. They have been obtained by taking the mf, — > 
limit of our rrib results written down in Sec. 8. For practical purposes the = 
results are sufficiently accurate for top decays since nib effects are generally quite 
small. This is particularly true if a running b quark mass at the top mass scale is used. 
Quantitative results on the mj, ^ corrections are given at the end of Sec. 8 as well 
as in 12^ 



7 Complete 0{cXs) results for nib = 

We are now in the position to put together our mf, = results. We add together the Born 
term results from Sec. 4, the one-loop results from Sec. 5 and the rrify — > limit of the 
integrated tree graph results according to Sec. 8. The mass and infrared singular log terms 
cancel among the 0{as) one-loop and tree-graph contributions as they must according to 
the Lee-Nauenberg theorem and one remains with a ffnite result. We choose to present 
our results in terms of scaled rate functions defined by Fj := Fj/Fg {i = U + L, + L^, 
U, L, F, S , U^, L^, F^, S^, A^) with Fq = Tu+L(BoTn) given by (x = mw/mt) 

To = Tu+L{Born) = \Vtb\ 5 . 34 

8v2 7r 

The angular decay distribution reads 



=^{l{t'u + PcosepTup){l + cos^ 6) + (35) 
4 7r 8 



d cos 9pd cos 9d(f) 4 tt 8 

+ ^(fi + Pcos^pfip)sin2 + 

3 - 

+ jiTF + PcosepTpp) cos9 + 



3 - 

H -i=TipP sin 9p sin 29 cos + 

2v 2 

3 - 

+ —^Vj!^pP sin 9p sin 9 cos ( 
v2 



where P is the degree of polarization of the top quark. As mentioned before one recovers 
the angular decay distribution written down in when substituting the mb = Born 
term expressions from Table 1 in Eg . (|35|) . 
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The various reduced rates f j are given by 



^ ^ 1 , [ (l-.'-^)(5 + 9.r^-6.r^) 2(1 -.r^)^(l + 2.r^)^^ 

27r ''(l-x2)2(l+2x2)\ 2x2 3x2 



(l-x2)2(5+4x2) 4(l-x2)2(l+2a; 



ln(l -x^) ^ '-^ ln{x) ln(l-a;^) -A{l+x^) x (36) 



„2M„,„^ 4(1-0:2)2(1 + 20:2)^^ 

9 



X (l-2x^) ln(x) - ^^^-^^^^^^U,{x^) , 



p _ l-2a;2 a ^2 f (l-a;)2(15 + 2a;-5a;2-12a;3 + 2a;^) 

~ 1+2^^2^ (1-^2)2(1 + 2x2) j 2^^ ^ 

(l+4x2)7r2 (1-^2)2(1-4x2), , ( 1 - x2) (3 - x2) (1 + 4x2) , , 

+ - — ^ H ^In 1-x ^In 1+x + 

3x^ x2 x^ 

4(l-x2)2(l-2x2)^_ ^ 4(2+5x^-2x6)^., ,1 
x^ x^ J 



1 + 2x2 27r " (l-x2)2(l + 2x2) [ ^ /V ■ / ■ 3 

„2 o™4\i„/„\ o^i 



+ in(l-x2) - 4(5+7x^-2x^) ln(x) - 2(l-x)^ x 



x^ 

(5+7x2+4x3), , 2(l+x)2(5+7x2-4x=^) , , 
X ^ln(x)ln(l-x) + — ^ln(x)ln(l+x) + 

X X 

2(l-x)2(5+4x + 15x2 + 8x3)^ ^ ^ 2(l+x)2(5-4x+15x2-8x3) ^ ^ xl, ^ 
- ^ — ^Li2 X +^ ^Li2 -X W38 

X X 



^ 1 x2 / (l-x2)(5+47x2-4x^) 27r2 

^ ~ T+2^2 + 2^^^(l-x2)2(l + 2x2) j 2^^ 3" ^ 

^ (l+5x2+2x^) _ 3(l-x2)2 ^ 16(1+2x2) ln(x) - 2(l-x)2x 

x^ x^ 

2-x+6x2+x3 2(l+x)2(2+x+6x2-x3) , , 
X ln(l-x)ln(x) — 5 ^ ln(x) ln(l+x) 



2(l-x)2(4 + 3x+8x2+x3) 2(l+x)2(4-3x+8x2-x3^ 



Li2(x) - ^ . ^Li2(-x) (39) 



X^ X^ 
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^2(l-.y(l+2.^) ^ 2(l-x^)(l-9x^+2.^) , 

+ %{\-x^f Li2(x) + 8(l+3a;2-a;^) Li2(-a;) I, (40) 



^ 1+2x2 ^27r ^(l-x2)2(l+2x2)\ 2x2 3x2 



(1-X2)2(2-5X2), 2^ 2^w N 4(1-X2)2 

+ —. ^ln(l-x^) -4(l-x^)ln(x) ^ — r-^ ln(x) In(l-x^) + 

x^ 



^^V^Li2(x^)|, (41) 



^ _ -2x2 ^ ^2 r_(l-x)2(12-55x+6x2-x3) _ 107r2 

1 + 2^2 + 2^ ^(l-x2)2(l + 2x2)\ 3~^ 

,^ ON 2(l-x2)2(l + 2x2) , , 2(l-x2)(7+21x2 + 2x^), , 
X 2+x2 +^ ^In 1-x +^ ^In 1+x + 

X^ X'^ 

+ 8(l-x2)2Li2(x) - 8(ll + 3x2+x^)Li2(-x) I, (42) 



la x2 f 

r^p = + — CiT- —, -\ - ( 15 -22x + 105x2 -24x=^+4x^) x 

^ 1+2x2 27r (l-x2)2(l + 2x2)\ ^ ^ 

^ (l-x)2 ^ (l+24x2 + 10x^)7r2 3(l-x2)2 ^^ _ (l-x2)(17+53x2) ^ 

2x2 3x2 ^2 2;2 

, ,^ , 4(l-x2)2 4(2 + 22x2 + 11x4)^ , 1 
xln(l+x)-^ — ^Li2(x) + ^ ^ ^Li2(-x)^, (43) 

^ 2x^ "s^ f / 2x/. 2x 2(l+x2 + 2x4)7r2 

4(l-x)2(l+3x+2x2+2x=^) 4(l+x)2(l-3x+2x2-2x^) , , , , , 
X — --{ — ^ — ln(x) ln(l+x) + 

X X 
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-Li2(x) + ^ — ^Li2(-x) , (44) 



X X 



^ I x^ { (l-x)2(ll-6x-7x2) (l + 2x2)7r2 

\p = \-—Ci? < — — + h 

^ l + 2x^ 2% ^(l-x2)2(l + 2x2)\ 2x2 3x2 

(l-x2)2(2-5x2) , , (l-x2)(2-9x2 + x4) , , 4(l-x2)2 

iAj JLj JLj 

X Li2(x) + ^^y^ Li2(-x)|, (45) 



-x a,^ x^ [(l-a;)2(12-7x + 12a;2) 7r2 



y2(l + 2a;2) 2ti " (l-a;2)2(l + 2x2) [ 6v^ 

(5 + 19a;2 + 2x^) (l-x2)2(l + 5x2) , , (l-x2)(l + 30x2 + 21x^) 

X ^ + ^ —pz ln(l — x) + ^ — — 1= X 

X 2^2x3 2^2x3 

, , , 2v^(l-x2)2 V2(7+15x2 + 4x^)^ , A 
xln l+x +^^^ ^Lia x -^-^ ^Lis -x , 46 

X X 



^ _ X X2 f (l-x2)(l+2x2) 7r2 

~ 72(1 + 2x2)^2^ ''(l-x2)2(l + 2x2) j 672 



X 



3-5x2 + 6x4 l-x2 2 i + 5a;2 X 5-11x2 , ^ ^ 
X ^ - '-^ ^ln(l-x2) ^ — ^ln(x) + 

(l-x)2(3 + 7x + 6x2) (l+x)2(3-7x + 6x2) , 
^ m(x) ln(l— x) -jzz m(x) m(l+x) + 

(l-x)2(7+15x + 10x2)^_ ^ (l+x)2(7-15x + 10x2)^_ ^1 



As mentioned in the Introduction the results for the total rate {U + L) agree with the 
analytical results given in ^, ^, |^ and in |T^. The six (mass zero) diagonal structure 
functions L, F and , L^, have already been listed in [|1^. They had been checked 
against the corresponding numerical results given in |2D|. The results on the non- 

diagonal structure functions and are new. As concerns the unpolarized transverse 
structure functions explicit expressions for the two linear combinations T+ = ^(U + F) 
and T_ = ^{U — F) relevant for the interpretation of the CDF measurement |]l| have been 
given in 
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We have also included 0{as) results on the unpolarized and polarized scalar structure 
functions Ts and Tsp- They determine the rrif, = unpolarized and polarized decay of 
the top quark into a charged Higgs {t —>■ b + H^) as it occurs e.g. in the Two-Higgs- 
Doublet-Model [2HDM). This can be seen as follows. The scalar projection of the 
Standard Model (SM) left-chiral current structure determines the coupling of the 

SM Goldstone boson, i.e. ijntPR — rribPi). This would be the coupling structure 

of the charged Higgs in the 2HDM when the ratio of vacuum expectation values is taken 
to be one. It is then evident that, for rub = 0, the scalar structure functions Ts and 
Tsp describe the decay t — > 6 + H'^ in the 2HDM irrespective of the value of the ratio 
of vacuum expectation values. The unpolarized scalar structure function Ts has been 



checked against the result of The result on the polarized scalar structure function 
Tgp is new. 

Before turning to the numerical evaluation of the various contributions we would like 
to discuss the large rrit limit of the various helicity structure functions. As expected from 
the statements of the Goldstone boson equivalence theorem the longitudinal and scalar 
contributions L, L^, S and 5*^ dominate in this limit. In fact, setting a; = one finds 



f,.f,^l + |c.(5-l.=). (48) 
fi, = fsP = l + ^Cf + (49) 

That Tl = ts and f^p = tsp for nit ^ oo can be understood from the fact that the 
longitudinal and scalar polarisation vectors e'^(O) and e^{S) become equal to each other 
in this limit since the longitudinal polarisation vector then simplifies to e''(0) = /m\Y + 
0{my/ / qo). The same observation is also at the heart of the proof of the Goldstone 
boson equivalence theorem. As concerns the tree graph contribution the statement that 
e^(0) = q^/mw + Oimw/qo) is certainly not true for all of three-body phase space as 
e.g. close to the phase space point where the W^-boson is at rest. The contribution from 
this phase space region to the three-body rate, however, becomes negligibly small when 
rrit — > oo. 

We now turn to our numerical results. As numerical input values we take mt = 
175 GeV and mw = 80.419 GeV. The strong coupling constant is evolved from as{Mz) = 
0.1175 to as{mt) = 0.1070 using two-loop running. The results are presented such that the 
reduced Born term rates are factored out from the reduced rates. This way of presenting 
the results allows one to quickly assess the size of the radiative corrections. One has 
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^U+L — 


1 - 0.0854, 


(50.a 




0.297(1 -0.0624), 


(50.b 


Tl = 


0.703(1 - 0.0951), 


(50.C 




-0.297(1 -0.0687), 


(50.d 


r'l'rr-Lr'iP 


0.406(1 -0.1162), 


(50.e 


' ^ = 


-0.297(1 -0.0689), 


(50.f 




0.703(1 - 0.0962), 


(50.g 




0.297(1 -0.0639), 


(50.h 


f,P = 


-0.228(1 -0.0810), 


(50.i, 


f^P = 


0.228(1 -0.0820), 


(50.j, 


= 


0.703(1 -0.0895), 


(50.k 


Tgp = 


0.703(1 -0.0922). 


(50.1 



The radiative corrections to the unpolarized and polarized rate functions are sizeable. 
They range from —6.2 % for Tu to —11.6 % for ^(_^^_^_J^^p compared to the rate correction of 
—8.5 %. The radiative corrections to the longitudinal and scalar contributions are largest. 
The radiative corrections all tend to go in the same direction. This is an indication that 
the bulk of the radiative corrections come from phase space regions close to the IR/M 
singular region where the radiative corrections are universal. When normalizing the rate 
functions to the total rate, as is appropriate for the definition of polarization observables, 
the size of the radiative corrections to the polarization observables is much reduced. 
For example, the 0{as) radiative corrections decrease the ratio Tl/Tu^l by 1.1% and 
increase the ratio Tu/Tu+l and the magnitude of the ratio Tp/Tu+L by 2.5 % and 1.8 %, 
resp., relative to their Born term ratios. The relative ratio Tu/Tl is increased by 3.6%. 
The values of the radiative corrections to the polarization observables are, however, large 
enough that they must be included in a meaningful comparison of future high precision 
data with the theoretical predictions of the Standard Model. 

The combination (Tu + Vf)/'^ determines the decay of an unpolarized top quark into 
a right-handed 1^-boson. This combination vanishes at the Born term level for mf, = 
as Eqs. (39) and (41) show. Adding up the corresponding numerical values of the 
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0{as) contributions in Eq. (51) one finds that the right-handed VT-boson occurs only 
with 0.094% probabihty. The 7^ effect in the Born term alone already amounts to 
0.036% (see Table 1). 

Altogether the 0{as) and the Born term mi, 7^ corrections to the transverse-plus 
rate occur only at the sub-percent level. It is safe to say that, if top quark decays reveal 
a violation of the Standard Model (SM) {V — A) current structure that exceeds the 1% 
level, the violations must have a non-SM origin. In this context it is interesting to note 
that a possible (V + A) admixture to the SM t ^ b current is already severely bounded 
indirectly to below 5% by existing data on 6 — > s + 7 decays ||35|, |37| . 

The rate combination {Tu + ^f)/'^ is in fact not the only combination that vanishes at 
the Born term level for = 0. Considering the fact that one must have Xw — Xt = —1/2 
at the Born term level the only surviving Born term level rate expressions are TZZ, Too 
and f Iq as alluded to before in Sec. 4. The notation employed for the reduced rates 
follows the notation used in Eq. ( p.3.a]) . The remaining rate expressions vanish at the 
Born term level but become populated at 0{as). They are 

f ++ = i(f[7 + + f^P + Tpp) = 0.000 833, 
foo=i(fz.-fip) = 0.000 389, 

f+- = {Tjp + Tap) = -0.000 236, (51) 
f;; = i(fc/ + - - tpp) = 0.000 093, 

fit = ^(Tu - f ^ + Tup - tpp) = 0.000 120. 



As remarked on before the latter two reduced rates Tz^'^ and Ttt vanish at the Born term 
level even for rrih 7^ since the net helicity of these transitions \Xw — At| = 3/2 exceeds 
that of the b quark lA;,! = 1/2. 

The four reduced rates r++, ToO) ^++ positive definite quantities since 

they result from squares of helicity amplitudes. Contrary to these T^q is an interference 
contribution and thus can be negative as it in fact is. In Eq. ( p]) we have also included 
the numerical values for the above five structure function combinations resulting from the 
(tree graph) corrections. They are all very small at the sub per mille level. 

In Sec. 4 (Fig.3) we have shown a lego plot of the Born term two-fold angular decay 
distribution in cos^ and cos6p. In order to be able to exhibit the size of the as corrections 
we show in Fig. 6 a contour plot of the same two- fold angular decay distribution with and 
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without radiative corrections, again setting P = 1. The radiative corrections are not very 
large in the upper two quadrants and become largest in the lower left quadrant of the 
contour plot when both cos^^ and cosOp tend to one. 

Instead of analyzing the three-fold or two-fold angular decay distributions one can also 
consider single angle decay distributions. They are obtained by integrating over the two 
respective complementary decay angles. For the cos 9 distribution one obtains 

= -{fu + 2rL){l + a0Cose + Pecos''e), (52) 

a cos 8 

where 



We have added the analytical and numerical Born term results for the asymmetry pa- 
rameters in brackets using x'^ = 0.211. The 0{as) values for the asymmetry parameters 
are ao — —0.357 and = —0.641, i.e. the ag corrections raise the magnitude of by 
2.3% and lower the magnitude of Pe by 1.5%. In Fig. 7 we show the cos^ distribution 
both for the Born term case and the radiatively corrected case. There is a pronounced 
forward-backward asymmetry. In the forward direction the differential Born term rate 
drops to zero. As discussed before the 0{as) rate does not vanish in the forward direc- 
tion due to real gluon emission. However, the radiative corrections are so small that the 
nonvanishing of the 0{as) rate in the forward direction cannot be discerned at the scale 
of the plot. In absolute terms the radiative corrections arc largest for cos 9^0 because 
of the large size of the radiative corrections to the longitudinal rate F^. Note that ag is 
not the conventional forward-backward asymmetry parameter which is defined by 

ciF(0>g>|)-dr(|>g>vr) _3 (_ 3 _ \ 

The Os corrections raise apB by 1.7% in magnitude. 
For the cos9p distribution one obtains 

dY^ 1 

^-{fu+L){l + Pae,cosep), (56) 



d cos 9p 2 
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where 



f . = 0.406Y (57) 



The tts-corrections lower agp by 3.4 %. 
Finally, the (p distribution reads 



where 



3tT^ Tap f 3tT^ X \ 

ld> = — = ^ = 0.597. 59 

The cos dependent contribution from Tj^ has dropped out because of having integrated 
over the full range of cos 6. If desired the contribution of Tjp to the (p distribution can be 
retained if one integrates only over half the range of cos^^. The a^-corrections raise 7^ by 
the small amount of 0.32 %. In Fig.8 we show the distribution both for the Born term 
case and the radiatively corrected case setting P = 1. 



8 Complete 0{cXs) results for nib 7^ 

Differing from the presentation of our = results in Sec. 7 we shall present our 7^ 
results in a form where each of the separate contributions to the rate remains identified. 
In particular we do not explicitly cancel the IR terms coming from the one-loop and tree 
graph contributions. We thus write 



J.QCD -p. ^ , mf\Vtb\^GF j ^ 1 , 4 as m^lVtbl'^GFx'^ 

^] X! /'W.* '^i"') + X! P{m.,n),i T^{m,n) + '^W + ^'^{m,n),i S{rn,n) \ 

y n=-l,0 m,n n=0,l rn,n ) 

The first term in Eq. (|60| ) represents the Born term contribution which is given by 

r.(Born) = %^fl^v^B. (61) 
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where the Bi are the Born term rates hsted in Table 1. The Born term contribution 
rj(Born) also appears as a factor in the third term where it multiplies the soft gluon 
factor «S'(A). The index i runs over the various structure function labels i = U + L, 
U'' + L^, U, ?7^, L, L^, F, F^, 5, , 1^ and . 

The second term in Eq. (|60D represents the one-loop contribution which is obtained 
by folding the one-loop amplitude in Appendix C with the Born term amplitude and then 
doing the appropriate projection onto the various structure functions. The appropriate 
coefficient functions listed in Table 2. The coefficient functions multiply the 

as one-loop amplitudes Fr = , F^ , F^ , Ff, F^, ^3"^ which are listed in Appendix C. 
We label the one-loop amplitudes consecutively by the index r = 1,...,6. Note that 
Table 2 contains only the vector current coefficient functions Ki^r (t = 1, 2, 3). The axial 
vector coefficient functions labelled by r = 4, 5, 6 can be easily obtained from the vector 
current coefficient functions by the substitution 



KpA 



KpV 
^1 



— Kl?V 



(62) 



The third term in Eq. (|60D contains the result of integrating the soft gluon function 
/S.SGF in Eq- (0)- The result depends on the (small) IR regularisation parameter A = 
nig/mt as indicated in the argument of the soft gluon factor S'(A). The universal soft 
gluon factor S'(A) is obtained by explicit integration and reads 



S{A) = -^CF\{l-x'+y') 
An 



Wi 



2 Li2(l - wiw^) + Li2(l - K) - Lis 1 - — + 



22/'ln(w;i) 



•-xy K' 



w,, 



+ \ lAw.w,) + Hn^^w,) { In (^) + i}] + 2 { In _ 2} + 



(63) 



where as in we have used the abbreviations 
X 1 — + — 



Wi 



y 1 + x^ - y"^ + \/X' 
In the limit y —>■ one has 



X 1 — + — 

y 1 + x'^ - y^ - y/X' 



(64) 



S{A) = -^Cp\il-x') 

4 TT 



4 + ln2|/-21nA + (l + 21nA)ln( ) + 

3 ^ y ^ 
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i 


Kpv,i 
1 






U+L 


v^((l-3/)2-a;2)((l+j,)2+2a;2) 


imt Va3(1+j/) 


imt VA3(l+y) 




A(l-2a;2-3/2) 


imtA((l+2/)2-a;2)(l-j/) 


|mtA((l+2/)2-a;2)(l-2/) 


U 


2V\iil-y?-x^)x^ 








C/^ 


-2x^X 








L 


yX((l-y)2-a;2)(l+2/)2 


imt VA3(1+j/) 


imt v^(l+y) 




A(l-?;2) 


|mtA((l+2/)2-a;2)2(l-y) 


imtA((l+j/)2-x2)2(l-j/) 


F 


-2 Ax^ 








pP 


2VA((l-i/)2-a;2)x2 








S 


VA((1+2/)2-x2)(1-2/)2 


imt VA((l+i/)2-a;2)x 


imt VA((l+y)2-x-2)x 


S"" 


A(l-?/2) 


imtA(l+x2-y2){l+j/) 


^mtX(l-x^-y'){l+y) 


7^ 












— ^7= mt V^x 
4 V2 


— !-= mt X 

4 V2 



Tabic 2: Coefficient functions Ki^^ that determine the contributions of the ag vector 
current one-loop amplitudes to the different rates {x = mw/mt, y = rrib/mt, A ~ 
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- 2 1- (rr^) + HI - In (y^) + Li,(x^) 



+ lnx2[>. (65) 



In agreement with the Lee-Nauenberg theorem the logarithmic dependence on the IR 
regularisation parameter A can be seen to cancel between the loop and the soft gluon 
contributions for each of the ten structure functions. 

The fourth term in Eq. (|60D finally contains the result of integrating the finite piece in 
the tree graph contribution Eq. (^), again after having done the appropriate projections. 
The result is given in terms of a set of standard integrals 7^(n), 'R-{m,n), <S(n) and S(rn,n) 
which are listed in Appendix A. Appendix B gives the values of the coefficient functions 
P(n),i, P{m,n),i, <7(n),j and that multiply the standard set of integrals in the various 

helicity structure functions. In Table 3 we have listed the range of values of the parameters 
m and n that characterize the different types of tree graph integrals. 

At this point it is perhaps appropiate to offer an excuse to the potential user of our 
nib results that our results are presented in a multiply nested form to be collected 
from Eqs.(|0|, |l|, H H, Table 2 and Append ices A, B and C. Contrary to the nib = 



results where a closed form representation was possible a presentation of unnested closed 
form expressions for rrib 7^ would require an extraordinary amount of space because of 
the presence of many different log and dilog functions and products thereof. Codes of the 
relevant expressions can be obtained from the authors on request. 

When we have evaluated Eq. ( pll| ) numerically the IR factors proportional to In A 
in the one-loop and tree graph contributions were set to zero by hand. The numerical 
evaluation of the remaining part is quite stable numerically. In particular the limit rrib — ^ 
is numerically quite smooth. This is demonstrated in Fig. 9 where we plot the bottom 
mass dependence of the total rate. Note that the 0{as) rate shows less dependence on 
the bottom mass than the Born term rate. 

The quality of the nib = approximation has been discussed before at the Born level. 
For example, at the Born term level the total rate is decreased by 0.27 % when going from 
nib = to nib = 4.8 GeV. Using the 0{as) nib ^ results from this section one finds 
that the 7^ corrections to the total 0{as) rate reduce the rate by 0.16% compared 
to the Born term reduction of 0.27%, i.e. the nib 7^ corrections to the as-contribution 
alone tend to counteract the 7^ effect in the Born term in the total rate (see also 
Fig. 9). The rUb 7^ corrections from the Os-contributions alone are surprisingly large 
considering the fact that the factor multiplying the a^-corrections Cp otsj (27r) = 0.023 is a 
rather small number. This can be understood in part by noting that the cts-contributions 
contain terms proportional to (mf/m^) ln(m^/m^) = —0.026 which is not a very small 
number. A further discussion of 7^ effects for the tts-contributions can be found in 
Il23| . Noteworthy is a large 20% correction to the 0{a^ transverse-plus rate r+ due to 
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i 


P\n) ' ^ 


P\m,n) ' ^ 






U+L 


— 


(-2,-1)... (0,-1) 


— 


(0,0), (1,0) 




-1,0 


(-2,0). ..(1,0) 


0,1 


(0,0),(0,1)...(2,1) 


U 


— 


(-2,1). ..(2,1) 


— 


(0,2)... (3, 2) 




-1,0 


(-2, 2). ..(3,2) 


0,1 


(0,0),(0,3)...(4,3) 


L 


— 


(-2,1). ..(2,1) 


— 


(0,2)... (3, 2) 




-1,0 


(-2, 2). ..(3,2) 


0,1 


(0,0),(0,3)...(4,3) 


F 


-1,0 


(-2,0). ..(1,0) 


0,1 


(0,0),(0,1)...(2,1) 


pP 




(-2,1). ..(2,1) 




(0,2)... (3, 2) 


S 




(-2,-1)... (0,-1) 




(0,0), (1,0) 




-1,0 


(-2,0). ..(1,0) 


0,1 


(0,0),(0,1)...(2,1) 


IP 


-1,0 


(-2, 2). ..(2,2) 


0,1 


(0,0),(0,3)...(3,3) 






(-2,1). ..(1,1) 




(0,2). ..(2,2) 



Table 3: Range of values of powers m, n in the different basic tree graph integrals 

mb effects That the bottom quark mass effect is so large in r+ can be appreciated 
in part by looking at the differential distribution in Fig. 5. We emphasize again that the 
mass effect may have been overestimated due to using a fixed pole mass, rather than a 
running mass which is smaller at the top mass scale. 



9 Summary and conclusion 

We have obtained analytical expressions for the 0{as) radiative corrections to the three 
unpolarized and five polarized structure functions that govern the decay of a polarized 
top quark. Although bottom quark mass effects are quite small in top quark decays we 
have retained the full bottom mass dependence in our calculation. In the limit uih —>■ 
the analytical results considerably simplify leading to compact expressions for the eight 
structure functions which are listed in the main text. The full mass dependence of our 
analytical results is written down in Sec. 8 and in the Appendices A and B. These finite 
mass results will prove useful for the theoretical description of 6 ^ c bottom meson and 
bottom baryon decays (see e.g. PB[). 

For top quark decays the radiative corrections to the structure functions range from 
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—6.2% to —11.6% where the radiative corrections to the unpolarized longitudinal struc- 
ture functions and the polarized structure function T(^^_^_x,)p are largest. These cor- 
rections are to be compared with the correction to the total rate which is —8.5 %. The 
radiative corrections to the structure functions all go in the same directions indicating 
that the bulk of the radiative corrections derive from contributions close to the IR/M re- 
gion of phase space where the radiative corrections are universal. Nevertheless the span of 
values of the radiative corrections exceeds 5 % and must be taken into account in a future 
comparison with precision experiments. The radiative corrections to rate combinations 
that vanish at the Born term level have been found to be rather small. In particular, the 
Os-correction to the normalized rate of an unpolarized top into positive helicity VT-bosons 
amounts to only 0.1%. As discussed in Sec. 7, the minuteness of the ag-contribution to 
positive helicity VT-bosons is of relevance when discussing a possible (V + A)-admixture 
to the Standard Model current. 

We have also determined the 0{as) corrections to unpolarized and polarized qi q2 
scalar current transitions. For t ^ b transitions these scalar current transitions are 
relevant for top quark decays into a bottom quark and a charged Higgs as they occur in 
the two-Higgs doublet model. For b —>■ c transitions these transition matrix elements are 
needed e.g. for the description of the semi-inclusive decays of the 5-mesons and the 
into spin- zero Dg mesons ^8|, pQ] . 

In this paper we have only studied the first order QCD corrections to the structure 
functions in polarized top decays. For the total rate one obtains a correction of —8.5 %. 



Second order QCD corrections to the rate are expected to amount to —2.6 % while 
electroweak corrections are known to increase the rate by +1.7% 5T]. For a high pre- 
cision comparison of theory and experiment of the structure functions it would therefore 
be desirable to calculate the two- loop 0{al) and the electroweak one-loop corrections 
to the eight structure functions. While the two-loop QCD corrections to the structure 
functions are very difficult and are therefore not likely to be done in the next few years 
the calculation of the one-loop electroweak corrections to the eight structure functions 
is presently under way Finite width corrections will also have to be accounted for. 
They lower the total width by 1.56% Q and affect the different partial helicity rates 
by differing amounts 



We would like to conclude this paper with a speculative note concerning a possible top 
quark mass measurement from an angular decay analysis using the fact that the structure 
functions are top mass dependent. This suggestion is much in the spirit of the suggestion 
of Grunberg et al. who advocated a similar measurement of heavy quark masses in the 



context of e"'"e~-annihilations [0. Assume that the percentage measurement errors on a 
L/{U + L) and U/(U + L) measurement are 5l and 5u, respectively. The percentage error 
on the mass measurement will be denoted by 5, i.e. we write mt = ffit{\ ± 5) where frit 
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is some given central value of the top mass. From the dependence of the respective Born 
term ratios on the mass ratio x = myy/rrit (assuming that the M^-mass is fixed) one finds 
that the percentage error on the top mass measurement is given by 5 = 5^(1 + 2xq) / (4xo) 
and hj 6 = 6u{l + 2xl)/2, resp., where we write x^ = xl{1^26) with xq = mw/^t- If we 
take rrit — 175 GeV as central value [xq = 0.211) this would imply that an 1% error on 
the angular structure function measurement would allow one to determine the top quark 
mass with 1.7% and 0.7% accuracy, depending on whether the angular measurement was 
done on the longitudinal (L) or on the unpolarized-transverse [U] (or for that matter 
(F)) mode. Since the radiative corrections change the ratios Vl/Vjjj^l and Vu/Vuj^l by 
1.1% and 2.4%, respectively, it is clear that one has to use the full 0{as) results for 
the angular structure functions if such experimental accuracies can be reached. This is 
illustrated in Figs. 10 and 11 where we plot the top mass dependence of Tl/Tjj^l and 
Tjj/Tuj^L for the Born term case and the 0{as) case for = 0. Note that the 0{as) 
curves are horizontally displaced from the Born term curves by approximately 3 and 3.4 
GeV, resp., meaning that one would make the correponding mistakes in the top mass 
determination from a measurement of the angular structure functions if the Born curves 
were used instead of the radiatively corrected ones. The present TEVATRON RUN I 
uncertainties on the top mass are around 4 % which is anticipated to be improved to 
1.7% during the initial stages of TEVATRON RUN II. It remains to be seen whether a 
mass determination based on angular measurements as proposed here can compete with 
the conventional method using invariant mass reconstruction. 
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A Integrals 

In this Appendix we catalogue the basic set of tree graph integrals that are needed in our 
nib calculation and give their analytical results. 



A.l Basic integrals 

In the first step of the tree graph integration one integrates over the gluon energy ko. 
After having done the integration on the gluon energy it proves to be convenient to 
perform a shift in the W-energy go integration variable by introducing the variable z — 
1 + — 2qQ/mt. One then encounters the following set of integrals 

f ,:"'(-/- f (Iz 
rl(m,n) '■= / / =, /^(n) := / / =, (66. a) 

{l-xr-e[ 



'(m,n) 



>(n) 



1, . l-x^+^+V^M) , (66.b) 



dz , A-x^ + z + ^\{l,x\z). 



,{z-y^) ^A"(l,x2,^) \-x^ + z-^\{l,x\. 



where A(l, a;^, z) = 1 -\- x'^ + z"^ — 2x?z — 1x^ — Iz. The required range of values of the 
parameters m and n are listed in Table 3. The cut-off parameters e'^ and 63 are needed to 
account for the spurious singularities which are artificially introduced by partial fraction- 
ing the integrands. The spurious singularities cancel as they must when all contributions 
to a particular helicity structure function are summed. 
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In order to get rid of the square roots the final substitution z —: l+x^—x{r + l)/r is 
introduced. The variable r has to be integrated in the interval [1+ei, 77 — £2], where 



V 



{1 + - + V\) /2x and X = X{l,x'^,y^) (67) 



as before. The spurious cut-off parameters ei and 62 replace the above cut-off parameters 
e'l and 62 and cancel in all final expressions. 

In order to keep our results at a manageable length we introduce the following set of 
auxiliary functions 



£,:=,„(1^), £.:=,„(!lli^), (68.C) 

TVo := Li2(r/a;) +Li2(^) -2Li2(x), M := Li2(r7x) - Li2 (^) (68.d) 

:= - ln(,) .„(1 + X) + in ((^p^^) !■> {^^i^) + (68.e) 

_Li,(l)+Li.((5^)+Li,(i^), 

^rj^ ^ rj — x ' ^ r\ — X ' 

Af, := - ln{rj) ln(l - x) - In ( l!L±lKi^) in ( ~ ^ \ + (68.f) 
_ Li.(- i) + L„((5-lil£) + L„(- 1^). 



r\ — X ' ^ rj — X 



and 



/3+(n) := (x - 1)" + (x + If, := (x- 1)"- (68.g) 

/5(n):=i^-il±^. (68.h) 

In the following we list our analytical results for the various types of integrals that are 
needed in our calculation. 
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A. 2 Integrals of type R (m,n) 



^(1,-1) 



7e(_2,-i) = ^ + - ^ Ar^, 69.a 

2 1 — x"' 2 

^(-1,-1) = -V/^ - (1 + ^')^^ + (1 - (69.b 

7^(o,-l) =^(l + x2-y2)AV2_^2^^^_^^^^ ^g9_^ 

-^A^/^ + (1 + x') f i(l + - y^)AV2 - x\C, - , (69.d 

=^-^rr^' = 2^5, (69.e 

7^(o,o) 7^(,o) = - ^, 7e(2,o) = - ^, (eg.f 

=(r^i^ + i-x^ 2 J' ^(-^'^) = i-x2 2 ' ^''-^ 

^(0,1) 7^(,l) = -AV^ + (l + x^)^^^, (69.h 

7^(2.l) =-^yV/^-^(l + x^)AV2 + (l + 4x^ + x^)^^, (69.i 

7e(3,i) = -Ix'/' + 1(1 + x')il + x'- y')\'l' - (3 + x'){l + 3x2) AV2 ^ ^gg ^ 

2\/-, I o 2 I 4\'^2 — -^1 



+ (l + a;")(l + 8a;" + a;^)- 



2 



7? _ 1 A 1 A 1 /I 1 \ . , 

/<'(-2,2) 4a;(i_a;)4 4a;(l + a;)'i (l-x2)2 1^1/2 (1 - a;)2y ' ^ 

^(1.2) =^^£3-2^A-^/?_(2)/:„ (69.m 



7^ 



(2,2) 





Ax 


(1 


-xf 




Ax 


(1 


-xf 


Ax 



^6 



Ax 2x 
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7^(4,2) =^^^^C,-^^^±^C,-^^^.{^)C,+ [{1- (69.P) 
X Q(l + X + - y'^f + (6 + 17x' + 6x^)) - 2(1 + x')A. 

A. 3 Integrals of type R („) 

= -A"= - (1 + - s')^^ + A'/Mn (^iV (70.a) 



A.4 Integrals of type S 

5(0,0) = - x\C2 - C,) - (71.a) 

5(1,0) = |(1 + + y^)AV^ - (2 + x')x'^:^^ - y'^, (71.b) 

5(0,1) = A^, 5(1,1) = (1 + x'Wo - A^/'>Ci + 2(1 - x')C, - ((1 - - y') , (71.c) 

5(2,1) = (1 + + a:^)Aro - ^(3 + 3x^ + y^)X'/^C, + 3(1 - a;^)£5 + (71.d) 
-\{{l-xr-y') {{l-xf + A + Sx^ + y'), 

Sio,2)^-^i^^2-^^3), (7i.e) 



2x 

(i + xi a-xi 



5(1,2) = -^^4^AA2 + ^-^-^N^ + M, (71.f) 



5(2,2) = --^-^^2 + ^^^^^3 + 2(1 + a;2)AAi + X^'^ - x\L^ - C,) - y'C,, (71.g) 
+ i(9 + 13x' + y')X'l' - (6 + 5x2)x2^^^ - y' (4(1 + x^) + y') ^, 
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'5(0,3) 


4x 1 — X 


'5(1,3) 




'5(2,3) 




'5(3,3) 





— X 1 + X X X J 

(1 - x)/:3 - (1 + x)C, - + ^A}, (71.j) 

X X 



1(2(1 - xf + (1 - x)=^/:3 - (1 + xfC^ - + ^i^A + 4xAro|, (Tl.k) 

4x XX-' 

1(2(1 - x)^ + (1 - xfC, - (1 + xfC, - + + (71.1) 

4x XX 

+ 12x(l + x^)Uq + 8x(l - x^)C^ - AxX^/'^Ci - Ax ((1 - xf - y^)) }, 

5(4,3) = ^{2(1 - xY + (1 - xyCs - (1 + xYC - + + (71.m) 

+ 24x(l + 3x^ + x^)A/'o + 28x(l - x^)C5 - 2x(7 + 7x^ + y'^)\^'^Ci + 
- 2x ((1 - xf - y^)) (7 + 9x2) ^ 



A. 5 Integrals of type S („) 

£2 \l/2 

'5(0) = + £1 In ( '-) + (£2 - -Ci) Iny + (72.a) 

+ L,(„.)-L,(^)-2L.((!^), 

5<.. - ^{- f ^. K^;^) - (- y - ^H- H^) }. (^^-M 

Of all the many integrals listed in A.2-A.5 the total rate calculation done before 
in ^, ^, ^, 1^, ^ requires only the five basic integrals 7^(_2,-i), 7?.(_i _i), 7?.(o,_i) , 15(0,0) 
and (5(1^0) compared to the 33 basic integrals that are needed for the full calculation. This 
may serve as a measure of the additional labour that is incurred when one calculates the 
complete set of structure functions as done in this paper. 



B Coefficient functions p(^), p(m,n), cr(n) and cr(^rn,n) 

In this Appendix we list the values of the various coefficient functions p(n),i, P(m,n),i, <^{n),i 
and a(^m,n),i (i = U+L, U^+L^, U, U^, L, L^, F, F^, S, S^, and A^) that multiply the 
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basic set of integrals listed in Appendix A as spelled out in the rate expression Eq. 
The coefficient functions involve polynomials in and which we sort by increasing 
powers of y"^. For reasons of conciseness we drop the suffix i denoting the particular type 
of structure function in the following listing. The contributions are collected in terms of 
powers of y'^. 

B.l Total rate i = U + L 

P(-2,-i) = , Uci.aj 

(l-x^)(l + 2x^) + (4-3xV + 3y^ 
P(-i,-i) = , (73.b) 

(3 - 2x^) + 

m-i) = -2 ' (73.c) 

^ y\{l + 2x^) + y^) 

^(0,0) = -2 , (73.d) 

{l + 2x^) + y' 

(J{ifi) = 2 , (73.e) 



B.2 Polarized total rate i = U^-\-L^ 



P(-2,0) = 


2/2(1 -a:2)2((l- 2x2) 




(74.a) 


x2 




P(-1,0) = 


(1 - 4x2 + 5x^ - 2x^) 


+ (5 - 4x2 - 5x^)?/2 - 2(3 - x2)|/4 


(74.b) 




x2 ' 


P(0,0) = 


2(1 - x2 - 6x^) - 3?/2 


+ y' 


(74.C) 


x2 




P(1,0) = 


(7 - 6x2) _ 




(74.d) 


x"^ 




p(-l) = 


^V^((l-2x2)-y2) 




(74.e) 


P(0) = 


oA((l-2x2)-y2) 
x^ 




(74.f) 


C(0,0) = 


yA((l-2x2)-y2) 

X^ 




(74.g) 
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cr(o,i) 


^4^2(1 - 2^2) (1 - x^) + (7 - - 6x^)^2 - (9 + + 

~ 2 ' 

X 


(74.h) 


^(1,1) 


^ (3 - + 6a;4) - (2 - x^)y^ - 

~ 2 ' 


(74.i) 




Jl-2x')-y''^ 

X"^ 


(74.j) 




^VX{l-x' + y'){il-2x')-y') 


(74.k) 


(J(l) 


X{l-x' + y'){{l-2x')-y') 

2 ' 
X"^ 


(74.1) 



B.3 Longitudinal rate i = L 



P(-2,l) = 


1/2(1 +y2)(i_:,.2)3 


(75. a) 


2 ' 
X"^ 




(1-^2) ((1 - x^y + {6 + x^- 3x^)y^ + (5 + Sx^)^^) 


(75.b) 




P(0,1) = 


(5 - 2a;2 - 7x^ + 4a;6) + (12 - SSx^ + a;^)y2 + (7 + x^)^^ 


(75.C) 


2 ' 




(7 - 31x2 + Ax^) + (10 + a;2)7/2 + 


(75.d) 


2 ' 


P(2,l) = 




(75.e) 


0-(0,2) = 


^y2((l + 10a;2 - lla;4) + (1 + a;2)2^2-) 

^2 ' 


(75.f) 


0-(l,2) = 


^ (1 + 10^2 - llx^) + (3 - 4^2 + + 2(1 + x2)y4 

2 ' 

x^ 


(75.g) 


(^(2,2) = 


^2(1 -3x2) + (3 + 2x2)^2 + ^4 

2 ' 


(75.h) 


0-(3,2) = 


2 ' 
X^ 


(75.i) 
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B.4 Polarized longitudinal rate i = L 



P(-2,2) — 


^2(i_^2)(i_^2)4 


(76 fi) 

I f \J •Cli J 






(1 - - + (7 - 2a;2 + 3x^)y^ - 4(2 + a;^)^^) 


(76.b) 






^2x^(3 - x^){l - x^ - {7- 5x^ + 6xV + (7 - ISx^ + 3x^)y^ 


(76.c) 


P(l 2) — 


^ (5 + 10x2 + 13x'^ - Ax^) - (9 + 48x2 + llx'^)y^ + 2(2 + x^)^^ 

9 ' 

X^ 


(76.d) 


P(2 2) — 


(16 - 20x2 - 4x^) - (17 + 20x2)y2 + y^ 


(76.e) 


X^ 


P(3,2) = 


x2 ' 


(76.f) 


P(-l) = 


o^(l-?/') 

° 2 ' 
X^ 


(76.g) 


P(0) = 


.A(l-|/2) 

X2 


(7D.nj 


C(0,0) = 


yA(i-i/2) 

2 ' 
X"^ 


(7f> \) 

I 1 yj.i J 


(Jen Q~i = 


^(1 - x2) (4x2(1 - x2)2 + (7 - 10x2 + 7x^ - 4x^)^2) 

x2 


(76.i) 


-1- 

1 


^(1 - X2)((9 + 8x2 _ 5^4)^4 _ 2(1 _ ^2)^6) 

2 ' 
X"^ 






^(3 - 5x2 ^ 17^4 _ 15^6) ^ - 9x2 ^ ig^4 ^ 112.6)^2 

x2 

^ (19 + 14x2 + llx^)y4 _ 4(1 + a;2)^6 

2 ' 


(76.k) 


(7(2,3) = 


^ (5 - 14x2 - 7x^) + (4 + 21x2 + llx^)y2 _ (n + 7x^)y'^ + 2y^ 

2 ' 
X^ 


(76.1) 


^(3,3) = 


^(l + 3x2) + 5x2y2_y4 

2 ' 
X^ 


(76.m) 


^(4,3) = 


2 ' 
X^ 


(76.n) 


(7(0) = 


v/A(l-|/2)(l_^2+^2^ 

9 ' 
X"^ 


(76.o) 
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a(i) = 4 , (76.p) 

B.5 Unpolarized- transverse rate i = U 

p^-2,i)^-2y\l-xy, (77.a) 

2(l-x^)((l-x^f -(l-5xV-2/), (77.b) 

P(o,i) = 2((l-6x2 + 5x^)-3(5-xV-2y'), (77.C) 

= 2{{17-5x')+y'), (77.d) 

P(2,i) = 2, (77.e) 

<7(o,2) = 4|/2((5-4x2-x^) + 2y2), (77.f) 

<7(i,2) =-4((5-4x^-x^)-2(2 + xV), (77.g) 

(T(2,2) =-4((6 + 2x2)+y2), (77.h) 

(^(3,2) = 4, (77.i) 

B.6 Polarized unpolarized-transverse rate i = 

P(-2,2)= 2y\l-xY, (78.a) 

= -2(1 - - + 2(1 + SxV - 2y'), (78.b) 

P(o,2) = -4((1 - x2)^(3 + x') - 2(1 + 3x^)y2 - 2(5 - x')y'), (78.c) 

P(i,2) = 4((9 + 10x=' + 5x^) - (27 + 5x')y2 + y4)^ (78_^) 

P(2,2) = 2(10(1 -x2) + 3y2), (78.e) 

P(3,2) = 6, (78.f) 

=-16VA, (78.g) 

P(o) = 16A, (78.h) 
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^(0,0) 


= 8\/A, 






(78.i) 


^(0,3) 


= 4(1 - x^){4:x\l - + (1 + _ 5^^.4^^2 , 


- 2(4- 




(78.j) 


^(1,3) 


= 4((1 - a;^)(3 - ILx^) + (9 - 22a;2 - llx^)y^ 


-2(1 


-2xV), 


(78.k) 


(7(2,3) 


= -4((13 + llx^) - (15 + 7x^)^2 + 2y^), 






(78.1) 


(7(3,3) 


= -4((l-5x^)+y2), 






(78.m) 


(7(4,3) 


= -4, 






(78.n) 


(7(0) 


= 8v^(l-a;' + |/'), 






(78.0) 


(7(1) 








(78.p) 



B.7 Scalar rate i = S 

_ y\l + y^){l-x^) 

P(-2,-i) , (79.aj 

_ {l + y^m-x-)+?,y^) 

P(-i,-i) - , (79.bj 

P(o,-i) = -3 ^2 ' (^9.c) 

(7(0,0) = -2 , (79.d) 

(7(1,0) = 2&±^, (79.e) 



B.8 Polarized scalar rate i = 

_ y\l-y^){l-x^f 
P(-2,o) , (80.a) 

(l-y2)((l-a;2)2 + 2(3-a;V) 
P(-i,o) = ^5 , (80.b) 

(l-y^)(2(l + 5x^)-y^) 
P(o,o) = , (80.c) 
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Pirn 


^2 ' 


(SO.d) 






(80.e) 


P{0) 


8A(1-|/^) 


(SO.f) 


^(0,0) 


^2 ' 


(80.g) 


(7(0,1) 


^ (1 - y^){4x\l - x^) + (7 + 5x2)1/2 _ 2^4) 

~ ^2 ' 


(SO.h) 


^(1,1) 


(1 - y^)(3(l - x^) + y^) 

~ 2 ' 


(80.i) 


(7(2,1) 




(80.j) 


(7(0) 


yA(l-y2)((i_^2)+^2) 

2 ' 


(SO.k) 


(7(1) 


A(l-y2)((i-a;2)+^2) 


(80.1) 


B.9 


Forward-backward-asymmetric rate i = F 




P(-2,0) 


9 /I 9\ 9 

= -2y2(l-a;2)2, 


(81.a) 


P(-i,o) 


n//i 2\2 1 /( 2 2\ 

= 2((1 - x^)"' + 4xV), 


(81.b) 


P(o,o) 


= 2(4(2 + x^) -7y^), 


(81.c) 


P{1,0) 


= -2, 


(81. d) 




/ — 

= 16\/A, 


(81.e) 


P(0) 


= -16A, 


(81.f) 


(7(0,0) 


= -8v^, 


(81-g) 


C7(0,l) 


= -4(4x2(1 - x2) + (1 + 5x2)|/2 - 2y^), 


(81.h) 


^(1,1) 


= -4(3(1 + x2)- 1/2), 


(81.i) 


(7(2,1) 


= 4, 


(81-j) 
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(7(0) ^-8V\iil-x')+y^), (81.k) 
(7(1) = 8X{{l-x')+y^), (81.1) 



B.IO Polarized forward-backward-asymmetric rate i = 

P(-2,i)= 2y\l-x^)\ (82.a; 

P(_i,i) = -2(1 - x^)((l - - (1 - hx^)y^), (82.b; 

P(o,i) =-2((l-6x' + 5x^)-(ll + xV), (82.c; 

P(i,i) =-2((ll + x2)+57/2), (82.d; 

P(2,i) = 10, (82.e; 

<7(o,2) =-V(l-a;2)(5 + x^), (82.f 

(7(1,2) = 4((l-x^)(5 + x2)-2xV), (82.g; 

a(2,2) = 4(2a;2 + |/2), (82.h; 

(7(3,2) = -4, (82.i; 



B.ll Polarized longitudinal-transverse-interference rate i = 

V2 y\l-xy 

P{-2,2) = ^ , (83.a) 

V2 (1 - x^iil - x'f + (3 + hx')y^ + 2y^) 
P(-i,2) = — 2 X ' ^ ' 

^/2 (1 - ^2)2(5 + 3x2) _ _ 38^2 + 29x^)^2 + 8(1 + x^)^^ 
P(o,2) = — ^ ^ , (83.C) 

(1 + 50^2 - 3x^) - (21 + 23x2)y2 + lOy^ 
P(i,2) = ^ ^ , (83.dj 

^ {b + 7x^)-2y^ 

P(2,2) = ^ ^ , (83.e) 
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=-4x/2^, (83.f) 

P(o) = 4V2-, (83.g) 

(7(0,0) = 2^2^' ^^^-^^ 

- x^)(4x^(l - x'f - (1 - + + 2(1 + 2x^)y^ 

cr(o,3) = v2 , (8d.ij 

/-(5- 18x2 + 5x^ + 8x6) -2(5 + + 8x^)^/2 + 2(1 + 4a;2)7/4 

(J(i,3) = V2 , (83.J) 

■^(3,3) =-v^;. (83.1) 

,,,, . 2V2 ^(('--')^^') . (83.m) 

X 
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B.12 Polarized parity-asymmetric rate i = 

V2y\l-x''f , 

V2 {i-x'){{i-xr + ^y') 

= — ^ , (84.b) 

V2 4(l-x2) + (3-7xV 

P(o,i) = ^ , (84.cj 

.(M, ^ (84.d) 

X — 2x'2 

^(2,2) = V2 , (84.g) 
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C Loop integrals 



In this Appendix we list the rrib 7^ one-loop amplitude corrections to the process t — »• 
b + W~^. They are determined from the vertex correction Fig. lb and the appropriate 
wave function renormalization constants Z2. We present our results in terms of the three 



vector current amplitudes 



1,2,3) and the three axial vector current amplitudes 



A 



1,2,3) defined in Eq. ( |25.a| ) in Sec. 5. Using the abbreviations in Eq. ( [6^ ) with 



F; 



T = f^w^ has 



^1' 



V 




mf + ml — 



m: 



2LU1 



w, - 



2112(1 



w,J 



+ 



+ i In (^-^) In(wiw^) + In (^) In ( 



w,. 



Wi 



-ln( 



2 2 



nit 



2g2 



ln(^)-4 + lnK^^)( 



m: 



'VX {nit + 



2g2 



j:,v ^ m - rrih j ^ fVfit + 2mb ml - ml 
q'^ I ^ mt — mf, g"^ ^ 



4tt 
ml^/X 



mb + (jrit — mb){3mt + mb) 



mt — mb 
3^ = F^{mt, mb) = F^{mb, mt) 



m 



+ 



i5) 



^6) 



^7) 



As before the IR singularity is regularized by a small gluon mass mg. The axial 
vector amplitudes F/^ can be obtained from the vector amplitudes by the replacement 
mt — > —mt, i.e. one has Ff^ijrit) = F^ {—mt) {i = 1,2,3). Our one-loop amplitudes 
are linearly related to the one-loop amplitudes given in The two sets of one- loop 

amplitudes agree with each other after correcting for a typo in mentioned in Sec. 5. 
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Figure 1 



Leading order Born term contribution (a) and 0{as) contributions (b,c,d) to 
t^b+W+. 




Figure 2 



— * 

Definition of the polar angles 9 and 9p, and the azimuthal angle 0. P is the 
polarization vector of the top quark. 
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Figure 3 

Born term Lego plot of the two-fold angular decay distribution 
dr/dcos9dcos9p with P = 1. 



10 



50 85 90 95 100 105 110 

qo [GeV] 



Figure 4 

Differential Vl^-boson energy distribution dVu+L/dqQ for the total rate resulting 
from 0{as) gluon emission (m^ = 4.8 GeV). 
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Figure 5 



Differential VF-boson energy distribution dr+ / dqo for the partial rate into pos- 
itive helicity H^-bosons resulting from 0{as) gluon emission for nib — 4.8 GeV 
(solid line) and for nib — (dashed line). 
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Figure 6 

Contours of the decay distribution of a fully polarized (P = 1) top quark in 
the cos 9p - cos 9 plane for mb — 0. The full lines are the distribution including 
the 0[as) corrections. 
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Figure 7 

Charged lepton polar angular distribution in the W rest frame for nib = 
(Born term: full line; 0{as): dashed line). Also shown are average values of 
the decay distribution. 
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Figure 8 



Azimuthal distribution of normalized rate for rub = (Born term: full line; 
0{as): dashed line). Also shown are average values of the decay distribution. 
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Figure 9 

Bottom mass dependence of the total rate Tu+l 
(Born term: full line; 0{as): dashed line). 
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Figure 10 

Top mass dependence of the rate ratio Tl/Vu^l for m;, = 
(Born term: full line; 0{as): dashed line). 




Figure 11 

Top mass dependence of the rate ratio Vu/^u+l for = 
(Born term: full line; 0{as): dashed hne). 
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